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Abstract. The ARGO-YBJ experiment has been in stable data taking for 5 years at the 
YangBaJing Cosmic Ray Observatory (Tibet, P.R. China, 4300 m a.s.l., 606 g/cm^). With a 
duty-cycle greater than 86% the detector collected about 5 x lO*^^ events in a wide energy range, 
from few hundreds GeV up to about 10 PeV. A number of open problems in cosmic ray physics 
has been faced exploiting different analyses. In this paper we summarize the latest results in 
cosmic ray physics and in gamma-ray astronomy. 


1. Introduction 

Aiming to face the open problems in Galactic cosmic ray (hereafter CR) physics through a 
combined study of photon- and charged particle-induced extensive air showers (hereafter EAS) 
with the same detector, ARGO-YBJ has been in stable data taking for more than 5 years at the 
YangBaJing Gosmic Ray Observatory (Tibet, P.R. Ghina, 4300m a.s.l., 606 g/cm^, 90.5° East, 
30.1° North). The detector acted simultaneously as a wide aperture (~sr), continuosly-operated 
7 -ray telescope at sub-TeV - TeV photon energies and as a high resolution GR detector in the 
broad energy range between few TeV and 10 PeV. 

The combination of the high elevation of the site, the full coverage of the central carpet and 
the high granularity of the readout provides very important advantages in various aspects of 
CR physics. In fact, in addition to the decrease of the threshold energy down to few hundreds 
GeV, the location of the esperiment above 4000 m a.s.l. ensures that shower fluctuations are 
small (we are working in the shower maximum region) and that all nuclei produce showers with 
nearly the same electromagnetic size. The low energy threshold is crucial for the overposition 
with direct measurements carried out by balloon/satellite-born detectors allowing one to cross- 
calibrate the different energy scales. The nearly independence of the size on the mass of the 
primary determinates the size - energy relation to be better defined, allowing the study of the 
elemental composition around the knee in a very efficient and reliable way. 

In this paper the latest results obtained by ARGO-YBJ in CR physics and in gamma-ray 
astronomy will be briefly summarized. 

2. The ARGO-YBJ experiment 

ARGO-YBJ is a full coverage air shower detector constituted by a central carpet ~74x78 m^, 
made of a single layer of resistive plate chambers (RPCs) with ~93% of active area, enclosed by 
a guard ring partially instrumented (~20%) up to ~100xll0 m^. The apparatus has a modular 
structure, the basic data acquisition element being a cluster (5.7x7.6 m^), made of 12 RPCs 
(2.85x1.23 m^ each). Each chamber is read by 80 external strips of 6.75x61.80 cm^ (the spatial 


pixels), logically organized in 10 independent pads of 55.6x61.8 cm^ which represent the time 
pixels of the detector [1]. The readout of 18,360 pads and 146,880 strips is the experimental 
output of the detector. In addition, in order to extend the dynamical range up to PeV energies, 
each chamber is equipped with two large size pads (139x123 cm^) to collect the total charge 
developed by the particles hitting the detector [2]. The RPCs are operated in streamer mode 
by using a gas mixture (Ar 15%, Isobutane 10%, TetraFluoroEthane 75%) for high altitude 
operation [3]. The high voltage settled at 7.2 kV ensures an overall efficiency of about 96% [3]. 
The central carpet contains 130 clusters and the full detector is composed of 153 clusters for a 
total active surface of ~6,700 m^. The total instrumented area is ~11,000 m^. For each event 
the location and timing of every detected particle is recorded, allowing the reconstruction of the 
lateral distribution and the arrival direction. The trigger is based on a time correlation among 
the pad signals depending on their relative distance. In this way, all the shower events giving 
a number of fired pads ^pad ^ ^trig in the central carpet in a time window of 420 ns generate 
the trigger. The whole system has been in stable data taking from November 2007 to January 
2013, with the trigger condition Ntrig = 20 and a duty cycle >86%. The trigger rate is ~3.5 
kHz with a dead time of 4%. 

Because of the small pixel size, the detector is able to record events with a particle density 
exceeding 0.003 particles m“^, keeping good linearity up to a core density of about 15 particles 
m“^. This high granularity allows a complete and detailed three-dimensional reconstruction of 
the front of air showers at an energy threshold of a few hundreds GeV. Showers induced by high 
energy primaries (> 100 TeV) are also imaged by the charge readout of the large size pads which 
allows to study the structure of the shower core region up to particle densities of 10Vm2 [2j. 

Details on the analysis procedure (e.g., reconstruction algorithms, data selection, background 
evaluation, systematic errors) are discussed in 0 0 E]. The performance of the detector 
(angular resolution, pointing accuracy, energy scale calibration) and the operation stability are 
continuously monitored by observing the Moon shadow, i.e., the deficit of CRs detected in its 
direction Pi- The measured angular resolution is better than 0.5° for CR-induced showers 
with energy E > 5 TeV and the overall absolute pointing accuracy is ~0.1°. According to 
MonteCarlo (MC) simulations the angular resolution for y-induced events results smaller by 
30% - 40%. The absolute rigidity scale uncertainty of ARGO-YBJ is estimated at 10% level in 
the range 1-30 TeV/Z 0[8]. 

3. Cosmic Ray Anisotropy 

The CR arrival direction distribution and its anisotropy has been a long-standing problem ever 
since the 1930s. The study of the anisotropy is complementary to the study of their energy 
spectrum and elemental composition to understand CR origin and propagation. It is also a 
powerful tool to probe the structure of the magnetic fields through which CRs travel. 

The anisotropy in the CR arrival direction distribution hs been observed by different 
experiments with increasing sensitivity and details at different angular scales. Current 
experimental results show that the main features of the anisotropy are uniform in the energy 
range (10^^ - 10^^ eV), both with respect to amplitude (10“^ - 10“^) and phase ((0 - 4) hr). 
The existence of two distinct broad regions, one showing an excess of CRs (called “tail-in”), 
distributed around 40° to 90° in R.A., the other a deficit (the “loss cone”), distributed around 
150° to 240° in R.A., has been clearly observed (for a review see, for example, [9]). In the 
last years different experiments reported evidence of the existence of a medium angular scale 
anisotropy in the both hemispheres [9]. 

In Fig. [1] the ARGO-YBJ sky map of medium angular scale (order of 10°) anisotropy in 
galactic coordinates as obtained with 4.5 years data is shown (for details see [TOl E]). The 
map center points towards the galactic Anti-Center. The zenith angle cut {6 < 50°) selects 
the declination region 6 ~ -20°-i- 80°. In this zenith angle bin the angular resolution is nearly 




Figure 1. ARGO-YBJ sky-map in galactic 
coordinates. The color scale gives the 
statistical significance of the observation in 
standard deviations. The map center points 
towards the galactic Anti-Center. 


Figure 2. Light component (p-|-He) en¬ 
ergy spectrum of primary CRs measured by 
ARGO-YBJ compared with other experimen¬ 
tal results. The results obtained by the 
ARGO-YBJ/WFCTA system are shown by 
the filled red squares [23] . 


constant and the energy of selected events is in the TeV range. According to the simulation, 
the median energy of the isotropic CR proton flux is ~1.8 TeV (mode energy ~0.7 TeV). 
No gamma/hadron discrimination algorithms have been applied to the data. Therefore, the sky 
map is filled with all CRs possibly including photons, without any discrimination. 

In spite of the fact that the bulk of Supernova Remnants (SNRs), pulsars and other potential 
CR sources are in the Inner Galaxy surrounding the Galactic Centre, the excess of CR is observed 
in the opposite, Anti-Centre direction. As stressed in [T2|, the fact that the observed excesses 
are in the Northern and in the Southern Galactic hemisphere, favors the conclusion that the CR 
at TeV energies originate in sources whose directions span a large range of Galactic latitudes. 
The right side of the sky map shows different few-degree excesses not compatible with random 
fluctuations (the statistical significance is up to 7 s.d.). The observation of these structures is 
reported by ARGO-YBJ for the first time, as discussed in m- 

From the theoretical viewpoint, the anisotropy of CRs is important as a direct trace of 
potential sources. Nonetheless, models failed to explain the whole set of observations at different 
angular scales [9]. Unlike predictions from diffusion models, that foresee only a dipole, the CR 
arrival distribution in sidereal time was never found to be purely dipolar, what translates in 
needing at least two harmonics to describe the data (the observation of anisotropies down to 
~10° implies multipoles of order 18). So far, no theory of CRs in the Galaxy exists which is 
able to explain the origin of these different anisotropies leaving the standard model of CRs and 
that of the local Galactic Magnetic Field unchanged at the same time. The anisotropy problem 
is the most serious challenge to the standard model of the origin of galactic CRs from diffusive 
shock acceleration m- 

4. Measurement of the Cosmic Ray Energy spectrum 

There is a general consensus that Galactic CRs up to the all-particle knee originate in SNRs 
accelerated by the first order Fermi mechanism in shock waves. The theoretical modelling of this 
mechanism can reproduce, in principle, the measured spectra and composition of CRs. Recently 
AGILE and Fermi observed GeV photons from two young SNRs (W44 and IC443) showing the 
typical spectrum feature around 1 GeV (the so called “vr^ bump”, due to the decay of vr*^) related 
to hadronic interactions mils]. This important measurement however does not demonstrate 















the capability of SNRs to accelerate CRs up to the knee and above. 

In the standard picture, mainly based on the results of the KASCADE esperiment, the knee 
is attributed to the steepening of the p and He spectra m- However, a number of experiments 
reported evidence that the bending of the light component (p+He) is well below the PeV and 
the knee of the all-particle spectrum is due to heavier nuclei (see, for example, [IZinBlllg]. 

A large number of theoretical papers discussed the highest energies achievable in SNRs and the 
possibility that protons can be accelerated up to PeVs (for a recent review see [2D] and references 
therein). Therefore, the determination of the proton knee, as well as the measurement of the 
evolution of the heavier component across the knee, are the key components for understanding 
origin and acceleration mechanisms of Galactic CRs. 

A measurement of the CR primary energy spectrum (all-particle and light component) in the 
energy range few TeV - 10 PeV is under way with the ARGO-YBJ experiment. To cover this 
wide energy range different ’eyes’ have been used: 

• ’digital readout’, based on the strip multiplicity, i.e. the picture of the EAS provided by the 
strip/pad system, in the few TeV - 200 TeV energy range [21]; 

• ’analog readout’, based on the particle density in the shower core region, in the 100 TeV - 
10 PeV range [2^ : 

• ’hybrid measurement’, carried out by ARGO-YBJ and a wide field of view Cherenkov 
telescope, in the 100 TeV - PeV region [23] . 

In the following we summarize results obtained below 100 TeV with the digital readout and 
describe the hybrid measurement up to the PeV range. Preliminary results obtained with the 
ARGO-YBJ analog data are described in [22]. 

4.1. Light component (p+He) spectrum of Cosmic Rays 

As described in [21], requiring quasi-vertical showers {6 < 30°) and applying a selection criterion 
based on the particle density, a sample of events mainly induced by p and He nuclei, with shower 
core inside a fiducial area (with radius ~28 m), has been selected. The contamination by heavier 
nuclei is found negligible. An unfolding technique based on the Bayesian approach has been 
applied to the strip multiplicity distribution in order to obtain the differential energy spectrum 
of the light component. 

The spectrum measured by ARGO-YBJ is compared with other experimental results in Eig. [2] 
(blue inverted triangles). The ARGO-YBJ data agree remarkably well with the values obtained 
by adding up the p and He fluxes measured by CREAM both concerning the total intensities and 
the spectral index [23]. The value of the spectral index of the power-law fit to the ARGO-YBJ 
data is -2.61±0.04. With this analysis for the first time a ground-based measurement of the CR 
spectrum overlaps data obtained with direct methods for more than one energy decade. 

This measurement has been extended to higher energies exploiting an hybrid measurement 
with a prototype of the future Wide Field of view Cherenkov Telescope Array (WFCTA) of 
the LHAASO project [25]. The telescope, located at the south-east corner of the ARGO- 
YBJ detector, about 78.9 m away from the center of the RPC array, is equipped with 16x16 
photomultipliers (PMTs), has a FOV of 14°xl6° with a pixel size of approximately 1° [26] . 

From December 2010 to February 2012, in a total exposure time of 728,000 seconds, the 
ARGO-YBJ/WFCTA system collected and reconstructed 8218 events above 100 TeV according 
to the following selection criteria: (1) reconstructed shower core position located well inside the 
ARGO-YBJ central carpet, excluding an outer region 2 m large; (2) more than 1000 fired pads 
on the central carpet; (3) more than 6 fired pixels in the PMT matrix; (4) a space angle between 
the incident direction of the shower and the telescope main axis less than 6°. This selection 
guarantees that the Cherenkov images are fully contained in the FOV, an angular resolution 
better than 0.2° and a shower core position resolution less than 2 m. 
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Figure 3. Scatter plot of the parameters 
Pc and pl for showers induced by different 
nuclei. The primary masses have been 
simulated in the same fraction, assuming a - 
2.7 spectral index in the energy range 10 TeV 
- 10 PeV. 



Figure 4. Light (p+He) component energy 
spectrum of primary CRs measured by 
ARGO-YBJ with different analyses. The 
systematic uncertainty is shown by the 
shaded area and the statistical one by the 
error bars. The parametrizations given in [281 
[25] are shown for comparison. A Horandel- 
like spectrum with a modified knee at Zxl 
PeV is also shown (for details see [22] 1. 


According to the MC simulations, the largest number of particles N^ax recorded by a RPC 
in a given shower is a useful parameter to measure the particle density in the shower core region, 
i.e. within 3 m from the core position. For a given energy, in showers induced by heavy nuclei 
Nmax is smaller than in showers induced by light particles. Therefore, N^ax is a parameter 
useful to select different primary masses. In addition, N^ax is proportional to where E^gg 

is the shower primary energy reconstructed using the Cherenkov telescope. We can define a new 
parameter p^ = logio{Nmax) — 1-44 • logio{Erec/TeV) by removing the energy dependence [23] . 

The Cherenkov footprint of a shower can be described by the well-known Hillas parameters 
1271, i.e. by the width and the length of the image. Older showers which develop higher in the 
atmosphere, such as iron-induced events, have Cherenkov images more stretched, i.e. narrower 
and longer, with respect to younger events due to light particles which develop deeper. Therefore, 
the ratio between the length and the width (L/W) of the Cherenkov image is expected to be 
another good estimator of the primary elemental composition. 

Elongated images can be produced, not only by different nuclei, but also by showers with 
the core position far away from the telescope, or by energetic showers, due to the elongation 
of the cascade processes in the atmosphere. Simulations show that the ratio of L/W is nearly 
proportional to the shower impact parameters Rp, the distance between the telescope and the 
core position, which must be accurately measured. An accurate determination of the shower 
geometry is crucial for the energy measurement. In fact, the number of photoelectrons collected 
in the image recorded by the Cherenkov telescope Npg varies dramatically with the impact 
parameter Rp, because of the rapid falling off of the lateral distribution of the Cherenkov light. 
Only an accurate measurement of the shower impact parameters Rp, and a good reconstruction 
of the primary energy allow to disentangle different effects. A shower core position resolution 
better than 2 m and an angular resolution better than 0.2°, due to the high-granularity of 
the ARCO-YBJ full coverage carpet, allow to reconstruct the shower primary energy with a 
resolution of 25%, by using the total number of photoelectrons Npg. The uncertainty in absolute 
energy scale is estimated about 10%. 

Therefore, in order to select the different masses we can define another new parameter pc = 
L/W — 0.0091 • (Rp/lm) — 0.14■ logio{Erec/TeV) by removing both the effects due to the shower 






























































distance and to the energy. The values of these parameters for showers induced by different 
nuclei are shown in the Fig. [3j As can be seen from the hgure, a suitable selection in the 
- Pc space allows to pick out a light composition sample with high purity. In fact, by cutting 
off the concentrated heavy cluster in the lower-left region in the scatter plot, i.e. pc ^ -0.91 
and Pc ^ 1-3, the contamination of nuclei heavier than He is less than 5%. About 30% of H 
and He survives the selection criteria. In the sample of 8218 events recorded above 100 TeV by 
the hybrid system, 1392 showers are selected in the (p-|-He) sub-sample. The aperture of the 
ARGO-YBJ/WFCTA system is ~170 m^sr above 100 TeV and shrinks to ~50 m^sr after the 
selection of the (p-|-He) component. 

The light component energy spectrum measured by the ARGO-YBJ/WFGTA hybrid system 
is shown in the Fig. El by the filled red squares. A systematic uncertainty in the absolute flux of 
15% is shown by the shaded area. The error bars show the statistical errors only. The spectrum 
can be described by a power law with a spectral index of -2.63 ± 0.06. The absolute flux at 400 
TeV is (1.79±0.16)xl0“^^ GeV“^ m~^ sr“^ s“^. 

This result is consistent for what concern spectral index and absolute flux with the 
measurements carried out by ARGO-YBJ below 200 TeV and by CREAM. The flux difference is 
about 10% and can be explained with a difference in the experiments energy scale at level of 4%. 
The measurement of the west-ward displacement of the Moon shadow under the effect of the 
geomagnetic field, as a function of the event multiplicity, allowed to calibrate the relation shower 
size - primary energy, thus calibrating the absolute energy scale of the detector at 10 % level 
below 20 TeV [ 6 ]. Above this energy the overposition with CREAM provides a solid anchorage 
to the absolute energy scale at few percent level. 

In order to extend the measurement of the ARGO-YBJ/WECTA hybrid experiment to the 
PeV, the selection cuts in the Pl pc space have been modified as follows: events for which 
Pl ^ -1.25 and pc ^ 1-1 are rejected. The aperture is a factor 2.4 larger. The contamination 
increases and the purity of the p-|-He sample below 700 TeV reduces to 93%. At 1 PeV the 
contamination is less than 13% increasing to 44% around 6 PeV. About 72% of p-|-He events 
survive the selection criteria. 

The preliminary measurement of the light component obtained by ARGO-YBJ with different 
analyses is summarized in the Eig. 01 As can be seen, all different approaches show evidence 
of a knee-like structure in the (p-|-He) spectrum starting from about 650 TeV, well below the 
all-particle spectrum knee confirmed by ARGO-YBJ at ~3 PeV [22]. With respect to a single 
power-law with a spectral index -2.63 the deviation is observed at a level of about 6 s.d. . 

More sensitive measurements of the energy spectrum and elemental composition in the range 
from 10^^ to 10^® eV are required to observe the knees of different nuclei and to investigate in 
detail the end of the spectrum of Galactic GRs and the contradictory results among different 
experiments. 

5. The Fermi Cocoon in the Cygnus region 

The Gygnus X region is one of the most luminous region of the Northern 7 -ray sky and it is 
rich in potential GR accelerator sites, e.g. Wolf Rayet stars, OB associations and SNRs. The 
Fermi collaboration detected a complex extended source, in a position consistent with the source 
ARGO J2031-I-4157, attributed to the emission of freshly accelerated GRs interacting with gas 
and radiation, filling a bubble (a “cocoon”) carved by stellar winds and multiple supernovae 
shock waves [30] . 

The significance map around ARGO J2031-I-4157 as observed by ARGO-YBJ is shown in 
the left plot of Fig. (5] To improve the sensitivity to 7 -induced events, an optimizated selection 
based on the shower core position is applied m- A smoothing is applied with an energy- 
dependent point-spread function evaluated for point 7 -ray sources. For comparison, the known 
TeV sources and the Gygnus Gocoon are marked in the figure. The sizes of markers indicate the 




Figure 5. Left plot: The significance map of the ARGO J2031+4157 region observed by 
the ARGO-YBJ experiment. Right plot: Spectral energy distribution of the Cygnus Gocoon. 
Different markers stand for the spectra measured by different detectors. The dot-dashed line 
shows the best fit to the Fermi-LAT and ARGO-YBJ data using a simple power-law function. 
The red solid (dashed) line is obtained by a hadronic model with a proton cutoff energy at 150 
(40) TeV. For references and details see [3T]. 


68 % containment size of the extension (for details see |31ji. 

The energy spectrum of the Cygnus Cocoon measured by Fermi-LAT, ARGO-YBJ and 
Milagro is shown in the right plot of Fig. [5j The flux determined by ARGO-YBJ appears 
consistent with the extrapolation of the Fermi-LAT spectrum suggesting that the emission of 
ARGO J2031-I-4157 can be identified as the counterpart of Cygnus Cocoon at TeV energies. 

The combined spectrum of Fermi-LAT and ARGO-YBJ can be described by a differential 
power law (dot-dashed line) dN/dE = (3.46T0.33) x 10“® • photons cm“^ 

s“^ TeV~^, suggesting the same origin for both GeV and TeV extended gamma-ray emission. 
Only statistical errors are shown, the systematic errors on the flux are estimated to be less than 
30%. The upper limits of Fermi-LAT and ARGO-YBJ indicate the presence of a slope change 
or cutoff below ~1 GeV and above ~10 TeV, respectively. At a distance of 1.4 kpc, the observed 
angular extension of about 2 ° corresponds to more than 50 pc, making the Cygnus Cocoon the 
largest identified Galactic TeV source. The observation of this large TeV source, not removed 
in the Milagro analysis of the 7 -ray diffuse flux from the Cygnus region, can explain the excess 
observed by them in the region 65° < 1 <85°, | 6 | <5°. 

Such a large region can be related to different scenarios. As discussed in [30], the favored 
scenario to explain the emission in the Cygnus Cocoon is the injection of CRs via acceleration 
from the collective action of multiple shocks from supernovae and the winds of massive stars, 
which form the Cygnus superbubble. Such superbubbles have been long advocated as CR 
factories, therefore the Cygnus Cocoon could be the first evidence supporting such hypothesis. 

In order to test a possible hadronic origin of the 7 -ray emission through the tt^ decay, we 
considered inelastic collisions between accelerated protons and target gas. We assumed that the 
primary protons follow a power law with an index similar to the gamma spectrum and with 
an exponential cutoff at 150 TeV, as suggested in [3D] to describe CR acceleration by random 
stellar winds in the Cygnus superbubble. This energy is the maximum proton cutoff allowed by 
the ARGO-YBJ upper limit. The resulting spectrum is shown in the right plot of Fig. [Djby the 
solid red line m- It is worth noting that the Milagro data are not described by this model and 
can be reconciled only with a cutoff of about 40 TeV (dashed red line) non consistent with the 










ARGO-YBJ results. 


6. Conclusions 

The ARGO-YBJ experiment has been in stable data taking for more than 5 years at the 
YangBaJing Cosmic Ray Observatory. With a duty-cycle greater than 86% the detector 
collected about 5x10^^ events in a wide energy range, from few hundreds GeV up to about 
10 PeV, exploiting the performance of the RPCs to image the front of atmospheric showers with 
unprecedented resolution and detail. 

Since November 2007 to January 2013 ARGO-YBJ monitored the Northern sky at TeV 
photon energies with a cumulative sensitivity ranging from 0.24 to ~1 Crab units, depending 
on the declination. This sensitivity, the lowest obtained so far, is allowing a detailed study 
of different galactic and extragalactic sources, focusing in particular on flaring and extended 
emissions. 

The high granularity of the readout and the high altitude location of the experiment offer 
a unique opportunity for a detailed study of several characteristics of the hadronic component 
of the CR flux up to about 10 PeV, starting from an energy window marked by the transition 
from direct to indirect measurements. In addition, the analog readout of the detector provides 
a powerful tool to study, for the first time, the particle distribution in the shower core region. 
Final analysis with the full statistics of the analog data will give new inputs to the hadronic 
interaction models currently used to describe particle physics and CRs up to the highest energies. 
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